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Abstract: Cosmological analysis based on currently available observations are unable to 
rule out a sizeable coupling among the dark energy and dark matter fluids. We explore a 
i variety of coupled dark matter-dark energy models, which satisfy cosmic microwave back- 

er I ground constraints, in light of low redshift and near universe observations. We illustrate 

^ I the phenomenology of different classes of dark coupling models, paying particular attention 

' in distinguishing between effects that appear only on the expansion history and those that 

\^ • appear in the growth of structure. We find that while a broad class of dark coupling models 

■ are effectively models where general relativity (GR) is modified -and thus can be probed 

I . by a combination of tests for the expansion history and the growth of structure-, there is a 

^ ! class of dark coupling models where gravity is still GR, but the growth of perturbations is, 

^ I in principle modified. While this effect is small in the specific models we have considered, 

H ' one should bear in mind that an inconsistency between reconstructed expansion history 

■ ■ ■ ' and growth may not uniquely indicate deviations from GR. Our low redshift constraints 

arise from cosmic velocities, redshift space distortions and dark matter abundance in galaxy 
voids. We find that current data constrain the dimensionless coupling to be |^| < 0.2, but 
prospects from forthcoming data are for a significant improvement. Future, precise mea- 
surements of the Hubble constant, combined with high-precision constraints on the growth 
of structure, could provide the key to rule out dark coupling models which survive other 
tests. We shall exploit as well weak equivalence principle violation arguments, which have 
the potential to highly disfavour a broad family of coupled models. 
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1. Introduction 

Cosmological probes [1-7] indicate that the universe we observe today possesses a flat 
geometry and a mass energy density made of ~ 30% baryonic plus cold dark matter and 
70% dark energy, responsible for the late-time accelerated expansion. The most economical 
description of the cosmological measurements attributes the dark energy to a Cosmological 

Constant (CC) in Einstein's equations, representing an invariable vacuum energy density, 
with constant equation of state w = —1. However, from the quantum field theory approach, 
the predicted energy for the vacuum fluctuations is ~ 120 orders of magnitude larger than 
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the observed value. This situation is the so-called CC problem. In addition, there is 
no proposal which explains naturally why the matter and the vacuum energy densities 
give similar contributions to the universe's energy budget at this moment in the cosmic 
history. This is the so-called why now? problem, and a possible way to alleviate it is 
to assume a time varying, dynamical fluid. The quintessence option consists of a cosmic 
scalar quintessence field which changes with time and varies across space, and it is slowly 
approaching its ground state e.g., [8-13]. The quintessence equation of state is generally 
not constant through cosmic time. In principle, the quintessence field may couple to other 
fields. In practice, observations strongly constrain the couplings to ordinary matter [14]. 
In addition, due to the the dynamical nature of the quintessence field, any coupling to 
the baryons would lead to time variation of the constants of nature, which are being 
tightly constrained see e.g., [15] and references therein. However, interactions within the 
dark sectors, i.e. between dark matter and dark energy, are still allowed by observations. 
A non-zero coupling in the dark sector could affect significantly the expansion history 
of the universe and the density perturbation evolution, changing the growth history of 
cosmological structures, see Refs. [16-28]. A number of studies have been devoted to 
analyze the constraints from Cosmic Microwave Background (CMB), large scale structure 
(LSS), Supernovae la and Baryon Acoustic oscillations (BAO) on a variety of interacting 
models [24,26,28-35]. Forecasts from CMB experiments, as the on going Planck-*^ and future 
satellite missions e.g., [36] on coupled cosmologies have been recently addressed [37]. 

Coupled cosmologies, in order to satisfy CMB constraints, predict values for the cosmo- 
logical parameters today which may differ substantially from the parameters values within 
non-interacting cosmologies. In order to fit high-precision CMB data available today, cou- 
pled cosmologies can hide their effects at very low redshifts. Therefore, low redshift probes 
are highly complementary and thus powerful to constrain interacting dark sector models. 
In this paper we focus on near-universe, low-redshift constraints in a variety of coupled dark 
matter-dark energy models. We explore the phenomenology of coupled models and consider 
what type of low-redshift observations are most suitable to improve present constraints. 
We pay attention in distinguishing between effects that appear only on the expansion his- 
tory and that can thus be tested with observations such as BAO and Supernova la and 
those that appear in the growth of structure. In the spirit of Rcf. [38] we shall exploit 
galaxy velocities and weak equivalence principle violation arguments as additional probes 
to tighten interactions among the dark sectors. We present as well unexplored, powerful 
constraints arising from dark matter abundance in voids, as well as from future, precise 
measurements of the Hubble constant Hq. Throughout this paper we assume a flat spatial 
geometry. 

It was suggested that the dynamical equilibrium of collapsed structures such as galaxy 
cluster could advocate in favour of an interaction between dark matter and dark energy [39]. 
The idea is that the virial theorem is modified by the energy exchange between the dark 
sectors leading to a bias in the estimation of the virial masses of clusters when the usual 
virial conditions are employed. This provides an additional near universe probe of the dark 

^http://www.sciops.esa.int/index.php?project=PLANCK 
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coupling. By comparing weak lensing and X-ray mass-observables to the virial masses, the 
authors of [40,41] claim that available data suggest that a small coupling could be present 
involving the decay of dark energy into dark matter. Notice that they also warn the reader 
that the lack of knowledge about the errors arising from systematics could weaken their 
conclusions. 

In Section 2 we propose a classification of interacting cosmologies into two broad 
families {DEvel and DMvel) and, in each of these families, two classes, depending on the 
scaling of the coupling with the matter or dark energy densities. This section presents also 
the background and the linear perturbation evolution for the different cases and discuss 
instabilities. Section 3 focuses on the phenomenology of these models compared to a 
ACDM model and to uncoupled models with arbitrary equation of state parameter for 
dark energy. Section 4 summarizes current constraints for interacting models and future 
prospects, devoting separate subsections to near universe Hq measurements, skewness tests, 
dark matter velocities and void dark matter abundances. Finally, we draw our conclusions 
in Section 5. 

2. Two times two families of deirk coupling models 

In this section we present the background and linear perturbation theory results for dark 
coupling models in general and then in the context of specific forms of the coupling. The 
models that we consider are naturally divided by two features of the model: whether the 
assumed momentum transfer is zero in the dark energy rest frame versus the dark matter 
rest frame. In each of these cases the the energy transfer could be taken to be proportional 
to pde or to 

2.1 Dark coupling models and modified gravity 

At the level of the stress-energy tensor it is always possible to introduce an interaction 
between the fluids of the dark sector in the following way [42] : 

VM?^fdm). = Q- and V^T5,), = -Q.. (2.1) 

The 4-vector Q y governs the energy-momentum transfer between the dark components 
and T(^ui),^ and T^^^^^^ are the energy-momentum tensors for the dark matter and dark 
energy fluids, respectively. Different expressions for the form of Qjy, that arise from a 
variety of motivations, can be found in the literature. Here we attempt to classify the 
couplings in broad families, based in the different phenomenology -from the astrophysical 
and cosmological point of view- they display. 

We consider two families of four momentum-energy transfer Q^^. In the first family of 
models (DEvel), the momentum exchange Qi, is parallel to the dark energy four velocity^ 

Qu = Q«i'^'V« (DEvel) (2.2) 

^The scale factor in the denominator is introduced because -at least at linear order- the velocity u is 
oc a, see Eq. (2.11). 
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In the second family of models (DMvel), Qi, is parallel to the dark matter four veloc- 
ity [21-24,26,27,38] ui"^^: 

Qu = Qu^^'^^/a (DMvel). (2.3) 

DEvel models oc Pdm^i'*^^ include all quintessence coupled models, see for instance, 
Refs. [10,16,23,32,38,43-51]. In DEvel models, there is no momentum transfer to the 
dark energy frame, so that momentum must be conserved in the dark matter frame. This 
implies a fractional increase in the dark matter peculiar velocity equal and opposite to the 
fractional change in energy density due to the presence of a coupling. This effect can be 
interpreted as an extra source of acceleration for the dark matter fluid, that will appear 
clearly in the dark matter velocity perturbation equation, see §3. By contrast, in DMvel 
models both momentum and energy density are transferred from the dark matter system 
to the dark energy one, and therefore the dark matter peculiar velocity field does not have 
this apparent force. 

The extra force effect in DEvel models should not come as a surprise: many quintessence 
coupled field models that appear in the literature can be written as a scalar-tensor Brans- 
Dicke [52] gravity theory. Let us also recall that f{R) gravity theories correspond to gener- 
alized Brans Dickc (BD) theory with a BD parameter wbti = or = —3/2, see e.g. [53] and 
references therein. The equations in the Einstein frame contain a new scalar field which 
satisfies Eqs. (2.1), being the energy momentum exchange Qy proportional to its 4- velocity 
(see e.g. [54]). 

The assumption of such a scalar interaction (DEvel) in the dark sector makes the 
acceleration of visible and dark matter particles different, inducing a "fifth force" effect 
(only for the dark matter), that is, a violation of the equivalence principle, for which the 
Eotvos experiments constraints do not apply. As noted by Ref . [55] , a large-scale fifth force 
in the dark sector might have substantial effect as a mis-match in the relative distribution 
of baryons and dark matter. DEvel models are therefore effectively "modified gravity" 
models. It is well known that deviations from the simpler ACDM paradigm in the form of 
general relativity (GR) modifications can be constrained following two approaches: a) at 
the background evolution level as inconsistencies between the high-redshift and the low- 
redshift universe [56, 57] and h) at the growth of perturbation level: a modified gravity 
model with the same expansion history as ACDM model, has a different growth of the 
dark matter structures. This has been extensively discussed in Refs. [58-60]. In fact, for 
some specific cases, there is a third possible approach: c) using weak equivalence principle 
violation (WEPV) constraints. 

In the following, we will see that option a) also applies to DMvel models and we 
illustrate its potential in the context of those models in §4.1. This approach only probes 
the background evolution and thus cannot be used alone to distinguish modified gravity or 
dark coupling from a minimally coupled dark energy model with general, time-dependent 
equation of state parameter w{z). For what concerns approach h), we show in §3 that the 
growth of perturbations in DMvel and DEvel coupled models can differ from the growth 
in an uncoupled model with identical background history. In the case of DEvel models 
(which in practice are effectively modified gravity models) it is already very well known 
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that growth can provide the key to break the existing degeneracy at the background level 
among GR and modified gravity. However, for DMvel models, the result we obtain here 
is a counterexample to what is commonly accepted in the literature: DMvel models are 
not a modification of gravity, but their growth can in principle still differ from the growth 
in a GR dynamical dark energy model which possesses the same background history since 
Pdm{z) docs not behave like dust (i.e., oc (1 + z)^). See also [61]. Finally, approach is 
specific of DEvel models (i.e., it is the smoking gun of those models) and we will consider 
it quantitatively in §4.4. 

2.2 Background evolution 

In Eq.( 2.2) and ( 2.3), Q drives the energy exchange between dark matter and dark energy. 

Indeed, one can easily show that for DEvel and DMvel families, the evolution equations 
for the dark matter and dark energy background energy densities reduce to: 

Pdm + ^npdm = Q , (2.4) 

Pde + 3H/9de(l +W) = -Q. (2.5) 

Pdm(Pdc) denotes the dark matter (dark energy) energy density, the dot indicates derivative 
with respect to conformal time dr = dt/a, H = d/a and w = -Pde/Pde is the dark-energy 
equation of state (P denotes the pressure). We work with the Friedman-Robertson- Walker 
metric, assuming a flat universe and pressureless dark matter = Pdm/Pdm = 0. The 
sign of Q determines the direction of energy transfer. For positive Q, the energy flows from 
the dark energy system to dark matter one. For negative Q, the energy flow is reversed. 
Note that if Q < 0, decreases with time because dark matter is being transformed 
into dark energy. The presence of a coupling Q also changes the dark matter and dark 
energy redshift dependence acting as an extra contribution to their effective equation of 
state. Indeed, the effective background equation of state for the two fluids are [24, 27] 

Therefore e.g., a negative Q yields an effective equation of state for dark energy that is 
more negative than w = Pde/pde- A negative Q will also contribute as a positive pressure 
in the dark matter background equation. Note however that the deceleration parameter 
satisfies, regardless of the presence of non-zero dark coupling, 

« = -^ = ^(l + 3^^de) , (2.7) 

where fide is the time dependent relative dark energy density. Therefore, we would still 
require w < —1/3 to have a universe with accelerated expansion. 

Because of the unknown nature of the dark sector, to-date there is no prescription in 
fundamental theory for a physically-motivated model for the coupling between the dark 
matter and dark energy fluids. The interaction term Q is currently mostly chosen in a 
phenomenological way. For the models considered in the literature in which the dark 



-5- 



coupling depends linearly on the dark sector energy densities, for each of the two families 
DEvel and DMvel, we propose the definition of two sub-classes of models: 

Q = Spdm (class I) (2.8) 
Q = Spde (class II) (2.9) 

where S denotes the interaction rate. The former classification should be considered as a 
"basis set" for the study of coupled models. For example, the coupling Q = Sipdm + ^2pde 
is also present in the literature [22,30,33,40,41] and it does not strictly belongs to class I or 
class II. Wc trust however that our analysis can guide the reader to determine what would 
be the relevant test to apply depending on the contribution that dominates at the scale or 
redshift corresponding to the probes under investigation. To gain physical insight on the 
phenomenology of coupled models we concentrate here on the simple cases of "pure" class 
I or class II models in both DMvel or DEvel scenarios. 

In coupled quintessence models S is proportional to the time derivative of the quintessence 
field, see e.g. Refs. [10,16,23,32,38,43-51]. More phenomenologically motivated models 
have taken E to be proportional to the Hubble expansion rate, since the former cosmolog- 
ical parameter has the appropriate time dependence ^. Even if the S oc ^ choice [21-24], 
is much more easy to handle mathematically than the S oc i?o case [21, 27] (being Hq 
the present-day value of the Hubble parameter), we will consider both possibilities in this 
paper. 

It is important to keep in mind that the expansion history docs not depend on the 
choice DEvel or DMvel (but does depend on the form of Q). Unfortunately, an analytic 
form for the expansion history H{z) cannot be written for generic Q. The only case 
where it can be done is in the case of Q = ^Hpde,dm (a mathematically easy to handle 
coupling) where both the coupling ^ and the equation of state parameter w arc constant 
(see Eqs. (A.l) and (A. 3) in Appendix A and e.g., [24]). Even if DMvel and DEvel models 
provide the same background history, the perturbation evolution is dramatically different, 
as we shall see in the next subsection. Therefore, geometrical probes alone are unable to 
distinguish among the two of them -or distinguish dark coupling from a uncoupled dark 
energy model with a w{z)-, even if, as we will see in §3, the two models are fundamentally 
different and different from minimally coupled dark energy models. 

2.3 Linear perturbations 

In the Newtonian gauge, the perturbed FRW metric at linear order in scalar perturbations 
is given by: 

ds'^ = a^[-(l + 2*)dr^ + (1 - 2^)dxidx'] , (2.10) 

where a is the scale factor, r is the conformal time, is the comoving coordinate, and * 
and $ are the scalar metric perturbations in the Newtonian gauge. 
The four velocity of a fluid reads 

u, = ai-il + ^),Vi) , (2.11) 
^For more discussion on a possible physical interpretation of this choice, see e.g., [21,22]. 
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where Vi is the fluid component pecuhar velocity. Since baryons are not coupled to the 
dark energy fluid, the continuity and the Euler equations for the baryons after decoupling 
are equivalent to those in uncoupled cosmologies: 

Sh = -Oh + 3i ; (2.12) 

= -ue^ + , (2.13) 

where 5 = Sp/p is the fluid energy density perturbation and 9 = div'^ is the divergence of 
the fluid proper velocity v^. We have assumed that baryons behave as a barotropic fluid 
with dPh/dph = 6Ph/6ph = 0, and we will assume the same properties for dark matter. For 
the coupled dark matter and dark energy components, we obtain, at linear order: 

-^dm = -(^dm + Sdm + ^) (2-14) 

Pdm \Q ) 

Odm = -nOdra + (1 " 6)-^(^dc " ^?dm) + fe'* (2-15) 

Pdm 



<5de = -(1 + w;)(0de - 3$) - ^ f ^ - 5de + * 

Pdc V Q 



-3W (Csde - ^) 



(2.16) 



^de = -n(l- 3cLe - ^de + T^&UMe + kH~h^^ 



Pdel + w' 

(2.17) 

We use the 6 notation introduced in Ref. [23], where 6 = refers to DEvel models with 
Qi, oc u^^^ (2.3) and 6 = 1 refers to DMvel models with Q^, oc ui^™^ (2-2); c^^^ is the 
dark energy pressure perturbation sound speed in the rest frame of dark energy. For the 
derivation of Eq. (2.16), we have used the following relation: 

SPde = C^sdc^Pde + (cLe " Cadc)Pdc^ , (2.18) 

which is valid in both Newtonian and synchronous gauges. If not otherwise stated, in 
the following, we assume c^^g = 1 (as for quintessence)^. The dark energy density per- 
turbations will not cluster significantly if the sound speed for the dark energy c^^^ = 1, 
and therefore, in the following, we can safely neglect dark energy perturbations in the 
perturbation evolution (see also the discussion in Refs. [25,38]). 

Notice that the Euler equation for the dark matter fluid is only modifled in the flrst 
family of models (DEvel) considered here, where Qi^ oc uif'^'^ case, violating therefore the 
weak equivalence principle, as we shall see in §4.4. This feature of DEvel models was 
previously discussed in §2.1. 

The growth equations for dark matter and baryons can easily be derived from Eqs. (2.12), 
(2.13), (2.14) and (2.15) going to the Newtonian limit (i.e. for k <C Ti). The growth equa- 
tion for dark matter is however rather sensible to the type of coupling. 

^Notice that current data still allows for ~ 0. However, in Ref. [62], it is shown that the two limits 
Cs ~ and 1 could be discriminated by future experiments. 
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The simplest case is the one of DMvel Q oc pdm (class I) coupled models for which we 
recover the standard perturbation equations at linear order as well as the growth equation 
for the two matter fluids a = dm, b: 



where q is the deceleration parameter of Eq. (2.7) and ' = d/da. Notice that the evolution 
equations for baryons are always the standard ones because they are not affected by the 
coupling introduced in Eq. (2.1). For the dark matter perturbations in DMvel Q oc pdm 
models, the difference with non-interacting cosmologies arise exclusively due to the dif- 
ferent background evolution of the quantities % and [25,38]. The growth equation 
is not modified but the scaling with redshift of pdm is different from that of a conserved 
pressureless fluid. As a result these models are effectively indistinguishable from minimally 
coupled dark energy models with a w{z). 

For all the other cases, i.e. DMvel Q oc pde class II and DEvel class I &l II, the dark 
matter linear perturbation equations and the growth equations are both modified by the 
presence of a coupling term Q. In order to derive the growth equation, an explicit form for 
the interaction term Q must be assumed, as we shall see in §3. 

2.4 Instabilities 

The existence of non adiabatic, early time instabilities in coupled models is a well known 
phenomena [21-24,27,63,64]. In brief, the dark coupling term which appears through p^^ in 
the dark energy pressure perturbations of Eq. (2.18) is a source for early time instabilities 
at large scales in the dark sector. The latter could arise in coupled models when the dark 
coupling term dominate on the standard (uncoupled) non adiabatic contributions to the 
dark energy pressure perturbation. The instability is then rapidly transferred to the others 
fluids components and to the curvature perturbation related to them, producing a non 
viable cosmological scenario. 

Here we briefly review the instability issue and provide a general recipe to avoid in- 
stabilities in the perturbation evolution in the case of a constant equation of state w (u^^^ 
can still vary). This recipe will be relevant throughout the paper, where the range and 
sign for the dimensionless parameters describing the coupling have been chosen to ensure 
an instability-free perturbation evolution. 

The onset of non-adiabatic instabilities depends on the form of the dark coupling Q 
(class I or class II), on the dark energy equation of state w and on the Qi, 4- velocity 
dependence (DMvel or DEvel). We shall define a doom factor d as [24] which is defined 
independently of the explicit form of the coupling Q: 



If |d| > 1, the interaction among the two dark sectors drives the non-adiabatic contribution 
to the dark energy pressure perturbation. At large scales, or equivalently, at early times 




(2.19) 



d = 



Q 



(2.20) 



3^Pde(l +1f) ' 
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("H/fc S> 1), the leading contributions in Q, or equivalently in d, to the second order 
differential equation for 5^^ reads: 



where 6 = (6 = 1) applies to DEvel (DMvel) models and ' = d/da. The sign of the 
coefficient of 5'^^ in Eq. (2.21) is, for c^^^ > 0, equivalent to the sign of the doom factor 
d. A positive doom factor d > 1 can lead to large-scale instabilities. Therefore, a simple 
recipe to avoid non adiabatic instabilities is to consider as viable models only those in 
which the energy transfer among the dark sectors Q and the (constant) equation of state 
of the dark energy component w are such that the doom factor d, Eq. (2.20), is always 
negative. 

Having defined all the tools, we briefly review the instability issue for the coupled 
models that will be analyzed in the following sections ^. We have shown that the key 
parameter is the doom factor, which only depends on background quantities such as Q. 
This means that the stability of the perturbations strongly depends on whether the model is 
a class I or a class II coupled model. However, the stability of the perturbations shows only 
a mildly dependence on the DEvel-DMvel classification of the model (see the h dependence 
of Eq. (2.21)). 

For class II models, it is easy to verify that, for Q = ^Tipde with ^ constant and 
negative, and w > —1, the doom factor is negative. In the following, we shall consider 

Qu = mpdeut"/a and Q, = ^HpdeutVa (2.22) 

with negative dimensionless coupling ^ and w > —1. The two models above are, respec- 
tively, the DMvel class II and DEvel class II models which we will analyze next. 

It was shown in Ref. [24] that for Q = ^^Pdm; the doom factor is always positive 
assuming w > —1 . It will also be enhanced by the Pdm/Pdc factor, see Eq. (2.20). A not very 
attractive solution to solve the problem would be to consider a constant w < —1 [22-24]. 
Another avenue is to use a time dependent equation of state for dark energy. In Ref. [27] 
it was shown that the DMvel model 

= -aTpd^u^^"^^ , with w{a) =woa + Weil-a) , (2.23) 

is stable for a specific range of constant wq, Wg and V/Hq. This is the DMvel class I model 
that we shall analyze. 

In the case of DEvel class I models, we shall refer to the case of a coupled quintessence 
model: 

Qu = /3pdmV^</./Mp (2.24) 

where (j) is the dark energy field driving the interaction, Mp denotes the Plank mass {Mp = 
I/VSttGtv) and /3 will be considered to be a constant. Note that in this case u'u^'' oc Vycp/cf). 

^For a detailed treatment of the instability issue in each of the cases reviewed below, we refer the reader 
to Refs. [21-24,27,63,64] 
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The early time non adiabatic instabilities discussed here arc different from the adiabatic 
instabilities. These adiabatic instabilities may arise at late times in quintessence coupled 
models if the sound speed of the total fluid gets negative [65] , even if Ref. [64] has shown 
that the slow-rolling of the quintessence field may avoid such instabilities in some cases. 
In the following, the term instabilities will refer exclusively to non adiabatic instabilities. 

3. Distinguishing dark coupling models from cirbitrary uncoupled w{z) 
dark energy models 

At the background level, a dynamical, redshift-dependent, effective equation of state w{z), 
can be mimicked by the combination of a constant equation of state w plus a dark matter- 
dark energy coupling Q [24]. Assume that the true, underlying cosmology possesses an 
interacting dark energy fluid with constant, non phantom, equation of state w > —1. If 
the background is analyzed setting the coupling to zero, one would reconstruct an effective 
redshift dependent w{z): 

_ 1 (1 + z)dRH/dz - mH 

where Rh{z) is a function of the background of the coupled model, 

Rh{z) = , (3.2) 

Hq 

and is matter energy density (dark matter plus baryons) as estimated from the CMB 
and extrapolated to 2; = assuming conserved pressureless dust. An analogous rela- 
tion between the redshift-dependent equation of state and the luminosity distance and its 
derivative was firstly presented in Ref. [66]. 

In the following, wc study uncoupled models with w(z) versus interacting DMvcl and 
DEvel class II cosmologies with Q = ^Hpde- We investigate how to distinguish among them, 
even if the background evolution in these cosmologies is identical. For this particular type 
of coupling, H{z) is given by Eq. (A.l), where we consider ^ < to avoid early-time 
instabilities. 

In the class II DMvel and DEvel models we consider here, the energy-momentum 
transfer is completely negligible at the high redshift relevant to the CMB. Assuming a flat 
universe and perfect measurements of fl^^h'^, i^y^h^, and the angular diameter distance to 
the last scattering surface from CMB observations, the amplitude of ^ is degenerate with 
the physical energy density in dark matter today ^. Figure I illustrates the values of ^^^^h? 
necessary to fit WMAP 5 year angular diameter distance data, i.e. the first acoustic peak 
position [1], as a function of the dimcnsionlcss coupling constant ^, taking into account 
that at the redshift of recombination z*, 0,dm{z*) — (^^dm + ^de^€/(3'"^e// ))(-*- '^^ 
equivalently, that the dark matter do not evolve as ^^mi^ + z)^ in the coupled models under 
study. The values of fide and H are defined coherently with the flat universe assumption. 

®There is a corresponding degeneracy with Ho that wiU be further quantified in §4.1. 
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Figure I: Present-day value of the dark matter energy density of the universe for the class II 
DMvel and DEvel models characterized by Q = S.'Hpdc as a function of the parameter ^, necessary 
to fit WMAP 5 year angular diameter distance data [1] and the physical dark matter and baryon 
densities at decoupling. The red solid (blue dashed) curve assumes an equation of state of the dark 
energy component w = —0.9 {w = —0.7). 



The values of the cosmological parameters and Hq, obtained as a function of the 
coupling ^, are presented in a Table in the Appendix B. The degeneracy direction depends 
only slightly on the (constant) value of w as shown in Fig. I. In the following, we shall 
restrict ourselves to ^ > —0.6, to ensure positive matter density at z = 0. 

It may be puzzling that dark matter could become negative, however, restricting our- 
selves to the ^ range above, this would happen in the future and by no means we can avoid 
it (the couplings ^ < —0.6 are already ruled out by a combined analysis of several cosmo- 
logical data [24,67]). An even more troubling feature is that in underdense regions Pdmix) 
may become negative well before ^dm does it globally. This statement depends on the 
smoothing scale employed to define Pdm{x)- Therefore, it may happen that at some special 
point in the Universe, for any value of the coupling, at some point in the past the current 
description makes the local dark matter density negative. However, notice that this would 
happen when fluctuations become highly non-linear, a regime which can only be explored 
with N-body simulations. Rather than an inconsistency of the model this indicates that 
the model should be seen as an effective, coarse grained description of some more complex 
and fundamental model. 

At the level of perturbations, the phenomenology in these coupled models is simpli- 
fied in the case where c^^ is large, where dark energy perturbations would be negligible. 
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Moreover we will assume in this paper that the perturbations in the expansion rate 6H can 
be neglected, see Ref. [28] for more details. In this situation, the rate of energy transfer 
between dark matter and dark energy is nearly homogeneous. For an energy density trans- 
fer, for instance, apdnu the product Pdm^dm remains constant, while S^m <^dm/(l ~ Q^)- 
Therefore we expect very small changes to the lensing signal, oc ridm^dm- However, since 
(5dm grows, we do expect some enhancement to the growth of structure in these models. 

In a more quantitative approach, one can derive the growth equation for dark matter 
in the Newtonian limit. We obtain: 

S^r. = -B^-f + I {a^J-^ + nj^^ , (3.3) 



where 



B = 2-q+{2-b)^^ , (3.4) 



Pdi 



m 



-C(l-g-3«;) + eM— + 1 



(3.5) 



3 ^^dm Pdm L VPdE 

Recall that 6 = for DEvel models and 6 = 1 for DMvel models. For negative coupling ^, 
the Hubble friction term B is suppressed and the A contribution to the source term in en- 
hanced. This implies that the dark matter growth in coupled models satisfying Eqs. (2.22) 
will be larger than in uncoupled models. More generally, this feature is valid for any cou- 
pled model in which Q is directly proportional to the dark energy density and Q/pde is 
negative, see e.g. Ref. [25]. 

One could be tempted to interpret the change in the source term of Eq. (3.5) as 
the result of a fifth force for dark matter, the total attractive force between dark matter 
particles being driven by some effective gravitational constant Ges with Ggs /G = A. This 
interpretation would mean the violation of the weak equivalence principle. However, such 
an interpretation is incorrect in the case of Q^, oc ul^^\ as the dark matter Euler Eq. (2.15) 
is not modified by the presence of the dark coupling and the velocities for dark matter and 
baryons are identical. 

We analyze now the linear growth in the three possible cosmologies that lie along the 
one-dimensional degeneracy defined by Fig. I. These three possible cosmologies, which 
have identical background histories, are: a) an uncoupled, albeit dynamical dark energy 
cosmology, with a varying equation of state w{z) given by Eq. (3.1), with H{z) given by 
Eq. (A.l) and assuming that Jldm in Eq. (3.3) scales with a~^, b) the coupled Q oc pde 
DMvel model with constant w = —0.9 given by the first equation of Eqs. (2.22) and c) 
the coupled class II Q oc pdc DEvel model with constant w = —0.9 given by the second 
expression of Eqs. (2.22). While in the DMvel model the baryon densities and velocities 
nearly trace those of the dark matter, this is not the case for the DEvel model. This 
difference will be carefully explored in §4.3-4.4. We consider observables probed by lensing 
{oc f^dm'^dm) if o^e ignores geometrical factors) and linear peculiar velocities (oc ^dm) see 
also § 4. 

Figure II shows the ratio of the values of f2dm'^dm Siud ^dm in cases b) and c) to their 
values in an uncoupled model with the same expansion history dictated by the relation in 
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Fig. I, for six possible values of the coupling ^. Let us emphasize that in Figure II, t he 
quantity f^dm<Jdm is derived using the solution for r^dm of Eq. (2.4) with non zero coupling, 
and ddm satisfying the growth equation (3.3). At the background level, the cosmological 
parameters were chosen to match the angular diameter distance to last scattering surface 
from CMB observations. i^dmjNC^dm.NC refers to the uncoupled model (Not Coupled) with 
the same expansion history than the coupled model. In the NC model H(z) is constructed 
assuming that ridm,NC has the usual dependence and the dark energy fluid is described 
by the effective equation of state of Eq. (3.1). The NC matter overdensity field (5dm,NC 
satisfies equation (2.19). An equivalent approach has been followed to derive 0dm/^dm,NC- 

For the DMvel class II model (case b)), there are only tiny changes to the velocity 
and lensing signals. While there is a clear divergence of the dark matter velocity and 
overdensity field in the DEvel class II model (case c)), we must be careful in considering 
how to apply observational constraints to this model. This enhancement of the velocity 
(for negative) yields the fifth force effect which is only present in DEvel models (see also 
the (1 — b) term in Eq. (2.15)). This effect will be described in more detail in §4.3.1. 

In summary, in both the DMvel and the DEvel models the dark coupling can, in 
principle, be distinguished from a generic, uncoupled, dark energy described by w{z) with 
perturbations growing as in standard GR. Note however that in the DMvel model GR 
is not modified. Although extremely small in this example, this class of models displays 
a mis-match between the reconstructed expansion history and the growth of structures, 
representing an exception to the commonly accepted interpretation that such a mis-match 
would be a tell-tale sign of modification of GR. In other words, if the measured expansion 
history and growth of structure are in agreement, this observation can be used to constrain 
deviations from GR. The conversely however is not necessarily true: a mis-match between 
expansion and growth could indicate deviations from GR or a dark coupling with GR 
unchanged. 

4. Low Redshift Observational Probes of Dark Coupling 

Coupled cosmologies, in order to fit the high-precision CMB data, predict values for the 
present-day cosmological parameters that differ substantially from the parameters values 
within non-interacting cosmologies. In other words, coupled cosmologies, in order to fit the 
high redshift (high z) observations, yield different expansion histories and different (shifted) 
values for present-day cosmological parameters. We show how to detect the induced shifts 
in the present-day value of some cosmological key parameters, as the Hubble constant Hq. 
Growth of structure probes are also explored and shown to be a powerful tool: while current 
errors on local bulk flows are still large, forecasted constraints on peculiar velocities offer 
a promising avenue. Let us mention that all the numerical results below were obtained 
choosing the cosmological parameters (as, for instance, r2(o)dm/i^) accordingly to WMAP 
5 year angular diameter distance data [1] as explained in section 3 for the Class II models 
studied here. For Class I DMvel models these numbers arc provided in Ref. [26]. Also, 
changes in e.g. last scattering surface, matter-radiation equality, and growth history for 
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Figure II: The ratio of the amplitude of dark matter peculiar velocities, 0dm, and the lensing 
signal, ndm^dm in coupled models compared to models with standard GR growth with identical 
expansion histories (labelled with NC for "Not Coupled"). The curves show ^ = 0,-0.1, ... — 0.6. 
The upper left panel shows the velocity ratio in the DMvel model and the upper right panel shows 
the lensing ratio. The two bottom panels show the same but for the DEvel model. 

cosmological perturbations are fully taken into account by means of a modified version of 
CAMB [68] which includes coupling effect at the background and linear perturbation level. 

4.1 Background High-z vs Low-z quantities mis-match: a worked example with 
Hubble constant measurements 

We first consider future CMB data (Planck) and future constraints on Hq (constraints 
from present cosmological data were already considered in Ref. [24] using the cosmomc 
package [69]). As it is often the case for forecasting the errors on parameters achievable 
with future data we make use of the Fisher matrix approach. 

Let us focus on the DMvel class II dark coupled models given by Q = iUpde (i-e. the 
first equation of Eqs. (2.22)) and calculate the shifts induced by the (neglected) coupling ^ 
on the present-day value of the different cosmological parameters. We apply the technique 
of Ref. [70] to compute the expected shifts at zero redshift on the cosmological parameters 
due to the presence of non-zero coupling in the data and of setting coupling to zero in 
the analysis. The authors of [70] developed a method which exploits the Fisher informa- 
tion matrix formalism. We briefly summarize their method in the context of interacting 
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models. Let us assume that we have two competing models: the M' model (a ACDM uni- 
verse) versus the M model (the DMvel coupled cosmology considered here). The method 
assumes that the n' parameters of model M' are common to M, which has p = n — n' 
extra parameters in it. These extra parameters are fixed to fiducial values in M' . For 
the form of dark coupling we illustrate here (the application to other coupled models is 
straightforward), p = 1, being the extra parameter the dimensionless coupling ^ of the 
DMvel interacting model, which is set to = in the M' model (ACDM). If the true 
underlying model turns out to be the M model (coupled cosmology), in the M' model, the 
maximum of the expected likelihood will not be at the correct values of the parameters. 
The n' parameters shift from their true values to compensate the fact that, effectively, the 
coupling ^ additional parameter is being kept at an incorrect fiducial value (i.e. is set to 
zero). If this (incorrect) fiducial value differs by S(f)^ from its true value, the others are 
shifted by an amount [70] 

66^ = -{F'-XpGp^6(t>^ , (4.1) 

where a,f3 = l...n', F' is the Fisher matrix for the model M' (ACDM) and is a 
submatrix of the M model (interacting cosmology) Fisher matrix. The set of n' parameters 
we use here to describe the M' model are the current baryon and dark matter energy 
densities Ob/i^ and ^dmh^, the current value of the Hubble parameter h, the amplitude Ag 
of the primordial scalar spectrum, the scalar spectral index and the equation of state 
of the dark energy component w. We have computed the shift in the n' parameters for a 
CMB experiment with the specifications of Planck The CMB Fisher matrix has been 
computed following Ref. [71] using a modified version of CAMB which includes the effect 
of the coupling at the background and linear order perturbation level. The Fisher matrix 
in the M' cosmology contains six parameters: ^dmh^ ,^bh^ ■> Hq, ng, w and Ag, being As is 
the amplitude of the primordial spectrum and rig the spectral index. The Fisher matrix in 
the coupled cosmology will contain an extra parameter, the dimensionless coupling ^. 

If we were living in a coupled Universe and future CMB Planck data were to be 
(wrongly) fitted to a ACDM cosmology, the value of Hq would shift by dlnHo/d^ = 0.3. 
This is the result of our Fisher analysis. However if the CMB data were to be fitted to a 
model where w is also a free parameter, new degeneracies open up and dlnHo/d^ increases. 
Thus if the Hubble constant Hq is determined independently with an uncertainty below 
3%, small values of coupling could be ruled out. A number of experiments (HST, Spitzer, 
GAIA and JWST [72]) are expected to measure Hq with 2% uncertainty in the next 
decade. Consequently, coupled models could be highly disfavoured over the next decade by 
the combination of CMB Planck data with precise measurements of the Hubble constant 
Hq. The other cosmological parameters will also be shifted from their true values but by 
amounts expected to be smaller than their combined statistical and systematic errors. 

The shifts for H{z) (at z > 0), as one may expect, are maximal at z = 0, thus 
forthcoming BAO observations which constrain H{z) with ~ % precision but only at a z 
where the volume element per unit solid angle is near maximal, will not offer significant 
further improvement. 

www.rssd.esa.int /PLANCK 
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To conclude, in this section we have provided quantitative forecasts of how future 
constraints on Hq can help in constraining the coupling between the dark sectors. 



4.2 Effects on the skewness 

It is well known that the skewness of the density field is relatively insensitive to the ex- 
pansion history (and thus the cosmological background) but is very sensitive to the growth 
of structure: mildly non-linear gravity leaves a distinct signature on the three-point cor- 
relation function (and therefore on the skewness) and galaxy bias also alters the skewness 
of the galaxy density field compared with the skewness predicted for the dark matter 
distribution. Here we consider the skewness as a test for the effects of dark coupling, fol- 
lowing Refs. [20,73]. A large effect on the skewness is a clear indication that the growth 
of perturbations is modified with respect to the standard expectation. Skewness tests are 
also affected by the ambiguity of whether galaxies trace the baryon or the dark matter 
distribution. 

Skewness is defined as the third order moment of the galaxy distribution 

o _ < ^i-r > (4.2) 



^ (< 6{xy >)2 ' 

where S{x) is the density contrast at the point x. In an Einstein de Sitter cosmology the 
skewness is predicted to be 34/7. This value has been shown to be independent of the 
nature of the dark energy component [74] or the background cosmology [75, 76] and it is 
only very mildly dependent on the shape of the primordial power spectrum (given current 
constraints on the spectral index). Here, we study possible deviations of the skewness from 
its standard value in the context of coupled cosmologies. Following Ref. [74] we compute 
the perturbation evolution in second order perturbation theory in coupled cosmologies, and 
compute the skewness parameter as 

53 = 4 + 26 = y + ^Q6-l) , (4.3) 

where e is a function of the first and second-order time-dependent components of density 
contrast {t) and 6^2) (*) > reads 

5(1) (t) is the solution of the growth equation at first order in perturbation which corresponds 
to Eq. (2.19) for uncoupled or DMvel class I models and to Eq. (3.3) for DMvel & DEvel 
class II considered here. Note, as we illustrate in the following, that the evolution equation 
for the second order density perturbation S(2){t) for the matter fluids a = dm,b depends 
as well on (5(i) {t) . 

Let us first analyze the case of Q oc pdm (class I) models. Ref. [73] has shown that 
the skewness deviations in DEvel class I cosmologies are ~ 1% for values of the coupling 
parameter satisfying current CMB constraints. For uncoupled and DMvel class I interacting 



- 16 - 



models, the equation for the second order density perturbation evolution for the matter 
fluids a = dm, b reads 



^(2),dm - 'i> a ^2 [ a2 ■ a2 / ' V a2 ' a2 

(4.5) 

For the DMvel class I interacting model of Eq. (2.23), the deviations of the skewness for 
both the cold dark matter and baryon distributions from its standard value is smaller than 
1% for the WMAP best fit of Ref. [26] and thus hard to be measured, even with upcoming 
future galaxy surveys. Recall that for DMvel class I models the continuity and Euler 
equations are exactly the same as for non interacting cases (after neglecting dark energy 
perturbations), and the only change in the growth is only introduced via the modified 
background evolution of 7i and pdm [25,38]. 

For the DMvel and DEvel coupled class II models with Q = S,7ipde of Eqs. (2.22), the 
second order growth equation for the dark matter fluid in the Newtonian limit is given by: 



^'(2), dm 3 



( ^^(2), dm ^ \ f(2)^ ^ 5(l),dm<^(l),b 



(2),dm - ^ a ^ 2 ^ ^ 

' '(4.6) 

where A and B are given by Eqs. (3.5) and (3.4) respectively. Let us recall that in order to 
obtain the growth equations for Q = ^"Hpdc models wc have neglected the contributions of 
5q. The additional terms proportional to ^6, ^5' in Eqs. (3.3) and (4.6) yield large effects 
on the cold dark matter distribution, and we observe ^ 10% deviations in the dark matter 
skewness. For the baryon distribution, the effects are below 1%. 

As already pointed out in § 3, for these two DMvel and DEvel class II models, the 
cold dark matter perturbation evolution can differ significantly from the baryon evolution, 
i.e. the changes in the growth of perturbations at linear and second order in these two 
models are not just due to a different background evolution. Even if the skewness effects 
in the cold dark matter distribution for these interacting models are large, it is not clear 
if they can be measured by galaxy surveys. Remember that in ACDM baryons and dark 
matter interact gravitationally in the same way at large scales so that the problem does 
not appear. 

If observations of the skewness of the galaxy distribution probe the baryon distribution, 
then skewness tests are not useful in constraining the dark coupling. Indeed, the presence of 
the coupling term only give rise to a 1% effect which is small compared with any foreseeable 
error-bar. If skewness observations (maybe via weak lensing, although recall that weak 
lensing probes the combination (5dm^dm and not d^m or alone) can yield the dark 
matter skewness, then this approach may be useful given the expected 10% effect. 

4.3 Peculiar Velocities 

The velocity field offers a test of the growth of structure, which is complementary to 
galaxy clustering, as, for example, is less sensitive to non-linearities and bias. Below we 
consider constraints that can be obtained from present bulk flow data in the local universe 
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and forecasted constraints achievable from future galaxy redshift surveys. These forecasts 
should be considered complementary to those presented in Refs. [24,26,28-35]. Throughout 
this section we will fix the CMB observables angular diameter distance to the last scattering 
surface and the physical dark matter and baryon density at decoupling. 

4.3.1 Constraints from peculiar velocities: Local bulk flows 

Recently Watkins et al. [77] have reported the observation of anomalously large bulk flows 
on 100 Mpc scales. In a Gaussian window of radius 50/i~^ Mpc, they find a coherent 
bulk motion of 407±81 km/s in conflict with the ACDM expectation of ~ 200 km/s at the 
2(7 level. Reference [78] pointed out that these results, if confirmed, would favour models 
with a growth of perturbations larger than ACDM model predictions. 

We first clarify the relation between the density and velocity fields in a general context 
where coupling is allowed. Wc again adopt a CMB-centric view: the amplitude of fluctua- 
tions at the epoch of recombination are well-measured, and so in the analysis that follows, 
we consider those to be flxed. We denote the amplitude of the density perturbation 
for a particular k-mode. The solution of the linear perturbations evolution (Eqs. (2.12) - 
(2.17)), which wc refer to by ^?dm/b(^) and (5dm/b(-2^)) describe the redshift evolution of the 
amplitude of each independent k mode: 

4,dm/b(^) = ~ T (4.7) 

t>dm/b(.^CMBj 

^ic,dWb(-) = A^J^^^^I^ , (4.8) 

o'dm/bl^CMBj 

where we have emphasized that when coupling is present, the growth of baryon and matter 
density and velocity perturbations are not necessarily the same. In the uncoupled case in 
the Newtonian limit {k S> W) and assuming that c^^j^ = 1, baryon and matter perturbations 
trace each other, and the velocity and density perturbations are related by: 

^k,dm/b = -a-H^ 4^ 7^r^ (uncoupled) . (4.9) 

'^dm/b(2^CMB) dhia 

In general, the relation between the growth equation solutions 5dm,b and Odm.b is not so 
simple. Nevertheless, one can eliminate the original mode amplitude ylk by defining 



? / ^ _ ^dm,b(^)/^dm,b(^CMB) , 
/dm,b(^) = -J —JJ -• (4.10) 

a-f2')dm,b(^)/ddm,b(^CMB) 

In the non-coupled case, / = d\n5 / d\n.a (the same for baryons and dark matter) , while 
in class II models {Q oc pdc) for both DEvel and DMvel, /dm = dhidfi^J dlna + ^p^^/ p^^. 
Finally, we can write the linear velocity field in terms of the density field as 

^ ifdm,haHSdm,h{k,a)k 
Vdui,h{^,a) = ■ p . (4.11) 
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Following Ref . [78] , the expected mean square velocity is 

{^^) = ^J^ dkk''P^{k)\W{k)\^ 

= ^ {A^Ai)\W{k)\\ 

27r^6'^(2:CMB) -'o 

where Pv{k) is the peculiar velocity power spectrum and W{k) is the Fourier transformed 
spherically symmetric window function In the second equality, the term outside the 
integral depends only on the normalized growth of velocity perturbations, 6{z)/9{zcmb), 
while the value of the integral depends on the power spectrum of initial fluctuations.^ This 
formula shows that bulk flows would be larger if the growth of perturbations is enhanced 
by some mechanism. In coupled dark matter-dark energy cosmologies, the growth factor 
can be modified in several ways [38] and therefore the predictions for the bulk flow may 
differ from those obtained in the context of a ACDM universe. Here we explore if coupled 
models could be favoured by bulk flow measurements. 

For the DMvel Q oc pdm interacting model of Eq. (2.23) both the continuity and Euler 
equations are unchanged and therefore, the expression for the uncoupled velocity amplitude 
/ = dlnS/dlna will still hold. 

Ref. [26] showed that in this model, CMB data alone disfavor positive values of F but 
still allow large negative values. For their CMB best fit model with T/Hq = —0.3, we 
obtain a bulk flow of 171 km/s. In principle, positive, higher values of F will predict larger 
values for the bulk flows. However, these values of F would not provide a good fit to CMB 
and/or BAO and SN data, see the analysis presented in Ref. [26]. Therefore, bulk flow 
observations do not indicate a preference for this particular coupled model. 

Figure III shows the bulk flow results for the DMvcl and DEvel class II models (Q oc 
Pdc) a-s a function of the coupling ^. The concordance ACDM cosmology predicts a bulk 
flow of (^2)1/2 = 203 km/s. As we have seen in S3, the DMvel model has dark matter 
peculiar velocities which are similar to those of non interacting models, since the Euler 
equation is unmodified. This is not the case for the DEvel model, in which the Euler 
equation changes, as can be seen from Fig. Ill: values of the coupling ^ < —0.35 lead to 
an effective 3 a level deviation from observations. 

However, this model must be handled with care since the dark matter and baryon 
perturbations are not equal. This large bias between the dark matter, baryon peculiar 
velocities and overdensities was shown in §3. Though smaller in amplitude, this feature 
is also seen in the coupled quintessence model explored with A?^-body simulations in Ref. 
[79,80]. While it is at first unclear which field the galaxies trace, we argue here that for 
^ > — 1 the velocities probed by, e.g., luminous red galaxies should be those of the dark 
matter, at least in the class II model considered in §3. It is generally accepted that the 

*For our numerical calculations, we use a Gaussian window function W{k) = e~*' ^ of radius R = 
50/1"^ Mpc. 

®This is not precisely true since (^k-Ak) is constrained tightly by the CMB in physical units {k in 
Mpc~^), while the radius of the window function is fixed at 50/i~^ Mpc. We account for this dependence 
on h in our calculations. 
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Figure III: The short (long) dashed curve shows the bulk flow predictions for the DEvel (DMvel) 
class II interacting model with Q = £,Hpdc- The cosmological parameters have been chosen to fit 
WMAP 5 year angular diameter distance data [1]. The dashed (dotted) blue lines represent the 
measured value of 407 ± 81 km/s, and the circle depicts the ACDM model prediction for the bulk 
flows (w2)i/2 ^ 203 km/s. 



massive luminous red galaxies (LRGs) were in place by redshift z > 1 (see Refs. [81-88] 
and references therein). At this redshift, the energy-momentum exchange is negligible. 
Therefore, galaxy formation should proceed as expected in ACDM - in the potential wells 
of dark matter halos. Moreover, we know that luminous red galaxies occupy massive dark 
matter halos today from weak lensing measurements [89]. 

To demonstrate that these galaxies remain bound in the potential wells of their dark 
matter halos, we need to demonstrate that the differential acceleration of the dark matter 
and baryons due to the WEPV "fifth force" in this model is much smaller than the force 
binding the galaxy to the dark matter halo. As discussed in §2.1, the choice of DEvel and 
Q oc pde implies a fractional increase in the dark matter peculiar velocity equal and opposite 
to the fractional change in energy density as energy is transferred from dark matter to dark 
energy. This corresponds to an acceleration on a dark matter halo with peculiar velocity 
V as 

dv\ 1 dm /A A\ 

— =v X — . (4.14) 



A dark matter halo has phalo ~ 200pcrit) so even for large couplings that correspond to 
depleting all of the mean dark matter density in the universe in the past 1 Gyr, 
1/m dm/dt < 1/200 Gyr~^. We take v to be the expected rms linear theory bulk flow, 
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~ 500 km/s, and find (|j)^gp^ < 8 x 10 m s ^. For a normal Newtonian orbit, the 
gravitational acceleration is 



dv\ GM _GMr^ir lOOi^^ryir 



I gravity '^eff '^vir ^ ^ 



(4.15) 



where M is the mass of the halo and refr is some fraction k < 1 of the virial radius ryir of the 
halo. For rvir ~ 0.5 Mpc appropriate for these galaxies, we get {^s) g-csNity > (^)wepv 
We conclude that under these assumptions, the galaxy peculiar velocities will trace the 
dark matter peculiar velocity field, and galaxy peculiar velocity observational results can 
be generally applied to interacting (DMvel or DEvel) class I and class II models. 

4.3.2 Constraints from peculiar velocities: redshift space distortions 

Even in the linear regime, peculiar velocities make the linear galaxy redshift-space power 
spectrum Ps(k) anisotropic when the underlying real space linear power spectrum i-*dm(^) 
is isotropic. These anisotropics go under the name of redshift-space distortions. In linear 
theory and under the flat-sky, distant observer approximation, we can generalize the re- 
sults of [59,90] to coupled models by simply replacing / = 9 In 5/9 In a with / defined in 
Eq. (4.10): 

<5f^(k, z) = (6gai + /V)<5dm(A:, z) , (4.16) 

where /x is the cosine of the angle between the wave vector k and the line of sight of a 
distant observer, 6gai the bias relating galaxy with dark matter overdensities in real space, 
i.e. 6gai = Sg/Sdim and 5dm{k, z) is the real space (isotropic) linear dark matter overdensity. 
Note that we assume in this section as well as in the previous one that galaxies trace 
the dark matter velocity fields rather than that of the baryons, since galaxies reside in 
dark matter halos. Therefore, in the notation of the previous section, galaxy redshift 
surveys constrain /(5dm = {^dm.{z) / 9 ^^{zc m b)0'H) {A\r^A1^^^'^ using the dependence in 
Eq. (4.16). Note that in the literature (e.g., [59]) the redshift space distortions constraint 
is often written instead as /erg, even though this probe is only sensitive to the amplitude 
of the velocity field, and not the amplitude of density perturbations directly. We follow 
this labelling convention, but replacing / — )■ / to account for the effects of dark coupling. 
Since coupled cosmologies can modify the growth of density and velocity perturbations in 
a significant way, one can use the measurements of redshift space distortions to constrain 
the strength of the dark sector interaction. 

We first consider the DMvcl class I coupled model of Eq. (2.23). For this form of 
coupling, as previously stated, both the continuity and Euler equations remain unchanged 
and in this case / = d\n.6/d\n.a. The only imprint of the dark sector interaction in the 
linear growth of perturbations for this first interacting model will therefore arise from the 
background quantities % and /9dm [25,38]. 

We also consider here the DMvel and DEvel class II interacting cosmologies of Eqs. (2.22). 
For these two cosmologies, the linear growth of the matter and velocity overdensities have 
extra contributions and therefore, the changes in the redshift space distortions observable 
will not arise exclusively from background quantities. 
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The results of § 3 indicate that the impact of couphng on velocity perturbations will 
predominantly be through the change in background in DMvel models, while sufficiently 
large couplings cause exponential growth of velocity perturbations in the DEvel model. 
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Figure IV: The red crosses denote the fas values for the model of Eq. (2.23), for T/Hq = —0.3 
and cosmological parameters fixed to model 2 of Ref. [26] . The blue triangles (squares) depict the 
fas expectations for a DMvel (DEvel) class II model with Q = ^fipdc for ^ — —0.5. The points in 
the left panel show current measurements of fas- The middle and right panels show the expected 
fas errors for BOSS and EUCLID-like surveys, respectively, assuming a fiducial ACDM cosmology. 



Figure IV illustrates the current and future constraining power of this method, where 
the red crosses show the fas predictions for the DMvel class I coupled model of Eq. (2.23), 
for T/Hq = —0.3, where the cosmological parameters have been fixed to their WMAP 
5 year data best fit values (see model 2 of Ref. [26]). The blue triangles (squares) refer 
to the predictions of the DMvel (DEvel) class II model explored here, with ^ = —0.5 and 
w = —0.9 (the remaining cosmological parameters are chosen accordingly to WMAP 5 year 
angular diameter distance data [1] ) . The black circles in the first panel of Fig. IV represent 
current fas measurements [91], in the other two panels they show the fas forecasts for 
BOSS and EUCLID-like galaxy surveys, assuming a ACDM fiducial cosmology. In order to 
forecast the errors from these two BAO galaxy surveys, we have used the Fisher information 
matrix formalism, combining the results from these two experiments with those expected 
from the on-going CMB Planck experiment. The Planck CMB Fisher matrix contains 
seven parameters 

ridm^^, ^^b/l^, -f^O, "-S, W, A, n-run, (4.17) 

where w is the dark energy equation of state, As is the amplitude of the primordial spec- 
trum, rig is the spectral index and rii-un is the spectral index running. The galaxy survey 
Fisher matrix contains one additional parameter, the growth factor /, plus a galaxy bias 
which is redshift dependent. We have marginalized over the galaxy bias in each redshift 
bin. 

BOSS-like survey parameters are the following: redshift coverage < z < 0.7, vlsky = 
10000 deg2 and a mean galaxy density of 2.66 x lO""^. For the EUCLID-like survey we 
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assume a redshift coverage of 0.15 < z < 2, a sky area of /sky = 20000 dcg2 and a mean 
galaxy density of 1.56 x 10~^. We have verified that our errors are fully consistent with 
those obtained by means of the Fisher routine provided by Ref. [92] . 

Firstly, note from Fig. IV that, for the values of the parameters chosen in this figure, 
the fas values for the DMvel class I model of Eq. (2.23) are much lower than current 
estimations. As carefully explored in Ref. [25], negative values of T (corresponding to 
positive Q) show a suppression of structure growth due to a smaller amount of dark matter 
in the past. The quantity fas will get depleted then, not only due to a lower growth factor 
/, but also due to a lower as- However, current errors on fas measurements are large and 
do not rule out these models at a high significance level, and other cosmological parameter 
choices could differ in the fas predictions. Secondly, note that the predictions for the 
DMvel class II model are very close to those of a non interacting model, as expected from 
the results presented on dark matter velocities for this model in §3 and in §4.3.1. Finally, we 
can observe from Fig. IV that for the DEvel class II model, huge fas values are predicted, 
especially at low redshifts, when dark energy starts to be dominant and the effect of the 
energy-momentum transfer starts to be important in the dark matter differential velocity 
equation. Such an effect can not be obtained in the two DMvel class I and class II models 
explored here because the Euler equation for the dark matter velocity is not modified in 
these interacting cosmologies. In the DEvcl class II model, couplings ^ < —0.4 will be 
ruled out at more than 3a with existing fas data. 

Accurate measurements of fas from future BAO surveys as BOSS or EUCLID may 
rule out significantly the form of couplings explored above, see the central and right panels 
of Fig. IV. For future surveys such as BOSS and EUCLID, the forecasted errors on the 
coupling for both the DMvel and DEvel class II models with Q = ^T-Lpdc can be as small 
as = 0.02 and < 0.01 respectively. These forecasted errors should be considered as 
optimistic, since in their derivation, we have assumed a perfect knowledge of the remaining 
cosmological parameters. 

4.4 Weak equivalence principle violation (WEPV) constraints 

Violation of the equivalence principle in the dark sector is a general feature of DEvcl 
models. Regardless of the underlying physics responsible for this effect, it can modeled 
phenomenologically by attributing to dark matter particles a "fifth force" specified by its 
potential. 

4.4.1 Coupled scalar field as an example of DEvel Q oc pdm models 

Kesden and Kamionkowski (K&K in the following) [93, 94], analyzed the consequences 
of WEP violation for dark-matter on galactic scales, focusing on dark-matter dominated 
satellite galaxies orbiting much larger host galaxies. It was shown by running several N- 
body numerical simulations that in the presence of a dark matter fifth force a pronounced 
asymmetry appears in the leading compared to the trailing tidal stream of the satellite 
galaxy. In summary, they inferred that a difference among dark matter and baryonic ac- 
celerations larger than 10% is severely disfavoured. They found that the leading-to-trailing 
ratio exceeds 0.5 for all simulations without a dark-matter force and never exceeds 0.2 for 
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a dark-matter force with 4% the strength of gravity. From the number of simulations per- 
formed, one can conclude that a difference between dark matter and baryon accelerations 
larger than 4% corresponds to at least a 2 (t deviation in leading-to-trailing ratio. However 
a more quantitative statement cannot be made as the way their simulations span the pa- 
rameter space for satellite interactions may not reproduce the same probability distribution 
as real interactions. 

To make this statement more quantitative we review the case DEvel class I models 
of the form of Eq. (2.24) considered in e.g., [79,80,95-99]. In these models the effective 
gravitational potential between cold dark matter particles is of the Yukawa-type. In the 
limit of massless dark energy interacting scalar field (f), the strength of the gravitational 
interaction is corrected by a contribution and reads: 

Geo = Giv(l + 2P^) (4.18) 

where Gn denotes the Newton's constant and /3 parameterizes the dimensionless dark 
coupling of Eq. (2.24) as well as the strength of WEPV. It should be stressed that in these 
coupled quintessence models, the correction to the standard gravitational potential term 

in the Eulcr equation for dark matter velocity perturbations is due to a contribution 
oc —fSk'^Scj) where 6(j) is the scalar field perturbation. Indeed, it can be shown [20], that in 
these models pScp ~ — /3^^, giving rise to a non negligible modification of gravity. Using 
Eq. (4.18), we see that the 10% limit from K&K reduce to |/3| < 0.22 in our notation. 

The authors of Ref. [100] obtained, in the notation used here, a 95% bound of < 0.11 
using CMB observations. More recently, it has been argued [101,102] that larger differences 
among dark matter and baryon accelerations, with an extra force as large as gravity, could 
be allowed if a screening length of = 1 Mpc/h is assumed. This screening allow the 
authors of Ref. [101,102] to evade cosmological-scales constraints (CMB). The caveat to 
this argument is that CMB constraints apply to redshift z ~ 1100 and the parameter ^ 
may vary with redshift. Ref. [103] find that large-scale structure observations at z < 0.2 do 
not show any indication for a typical scale where this Yukawa-type modification switches 
on, but constrain the strength of the modification only weakly; [104] improve these bounds. 

A possibly promising avenue to be explored from forthcoming large, SZ-selected galaxy 
clusters surveys, is that of the baryon mass fraction in galaxy clusters. In fact, as shown 
in [80], a WEPV , with increasing coupling reduces the clusters baryon fraction: a /9 = 0.2 
reduces /b in clusters by 10% with respect of a ACDM model. This effect is mostly due to 
a mis-match of velocities on large-scales not driven by Eq. (4.18). Tantalizingly, this could 
alleviate the tension between cosmological and clusters determinations of /b from current 
data, however to-date the error-bars are still large. The results of [80] derive from N-Body 
simulation for coupled models. This type of analysis had already been first performed 
in [79] and [80] claims to obtain similar results on the enhancement of the bias in the 
nonlinear region within and around massive halos. On a another hand, let us mention 
that the results of [80] completely disagree with those of [79] on the increase or decrease 
of the inner halo overdensity with respect to CDM assuming the existence of a coupling 
between the scalar field and dark matter. More constraints have been derived on these 
type of coupled models using background up to cosmological perturbations evolution (see 
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e.g. [98,99]), and can affect the late-time Integrated Saclis- Wolfe (ISW) effect at large 
scales [105]. See also [106] for constraints on time dependent coupling. 

Let us repeat that on cosmological scales, a WEPV only in the dark sector, would 
yield a mis-match between the baryon (and thus galaxy) distribution and the dark matter 
one [102] which is scale-dependent (see also [80]). By comparing left (weak lensing) and 
right (galaxy clustering) panel for fig. 2 of [104] we conclude that there is no evidence 
for a large mis-match from present data, however the error-bars are still large. Forthcom- 
ing gravitational lensing surveys, combined with galaxy and Lyman alpha surveys should 
therefore constrain these models further. 

So far we have considered very specific models of the DEvel family for which a po- 
tential can be specified. The above consideration hold qualitatively also for the more 
phenomenologically-defined models as we show next. 

4.4.2 DEvel Q oc pde model 

As an example, we illustrate also the effect of from WEPV for the DEvel class II toy-model 
considered along this paper and characterized by Qy = ^Hp^^ut^- In this DEvel class 
II model still the dark matter and baryon accelerations are clearly different, i.e. compare 
Eq. (2.15) with 6 = to Eq. (2.13). The results obtained in the following for this particular 
choice of coupling can be easily generalized to other choices of interaction. At sub-horizon 
scales, in the Newtonian regime {k 3> H)-, the baryon and dark matter accelerations are 
given by: 

= HSb + A;2*. (4.19) 

^dm = n(l+ f (jdm + mSdm—] + fc'* , (4.20) 

where we have neglected S^e and SH contributions and A;^* reads 

k^^ = An\n,,6^ + J^dm^dm) • (4.21) 

Figure V shows the effect of the dark coupling on the relative dark matter-baryon accel- 
eration for the model illustrated here, for two different assumptions about the dark energy 
equation of state w. The cosmological parameters (as, for instance, r^dm^^) have been 
chosen accordingly to WMAP 5 year angular diameter distance data [1] as explained in 
section 3. 

There is clearly a well defined relation between the strength of coupling and the mis- 
match between the dark matter and baryon acceleration that is robust to changes in w and 
is well approximated by a power law: 

In ^^Lzi^^ = 2.69 + 0.977 ln(-0 (4.22) 

up to mis-matches of order of few x 100%. This has been derived assuming linear theory 
and therefore quantitatively strictly applies only to cosmological linear scales. 
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Figure V: Relative dark matter-baryon acceleration for Qi, — ^Hpdculf"^ /a using Eq. (4.19) 

and (4.20). The cosmological parameters have been chosen to fit WMAP 5 year data, see text for 
more details. The red solid (blue short dashed) curve assumes an equation of state of the dark 
energy component w = —0.9 {w = —0.7). 

From the discussion above, our numerical results seems to indicate that for DEvel 
class II models with = ^Tip^eut^ /a, couplings ^ smaller than ~ —0.002 would lead to 
a 2fT sigma deviation in the leading-to-trailing ratio and couplings smaller than ~ —0.01 
are strongly disfavoured by K&: K numerical analysis. Let us mention that those results 
are to be taken with care given that no asymmetry in leading compare to the trailing of 
e.g. Sagitarius dwarf galaxy is observed while K&: K numerical best fit model tends to a 
leading-to-trailing ratio of 0.66 in the absence of any dark matter fifth force 

4.5 Constraints from matter abundance in galaxy voids 

An uncontroversial (indisputable) feature of these models arises from the fact dark energy 
is smoothly, uniformly distributed but dark matter is not, and most of the volume of 
the Universe is occupied by regions under-dense of dark matter (voids). The coupling 
Q oc pde characteristic of class II models requires that, independently on the local dark 
matter density, a given amount of dark matter per unit volume must turn into dark energy. 
Thus in the underdense regions the coupling Q oc pde must be just an effective, coarse- 
grained description: in fact, once locally in an underdense region all dark matter has been 
transformed into dark energy, the description adopted must break down. In addition, the 
uniform depletion of dark matter makes the overdensity in voids approach 5 ~ — 1 very 
rapidly, implying a prompt breakdown of linear perturbation theory. While a detailed 
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modeling of this behavior may require numerical simulations, it is clear that observations 
of underdense regions offer a promising avenue to constrain such models. We consider now 
observations of voids properties to constrain these models. 

Coupled dark energy-dark matter models lead to a dark matter depletion (or en- 
hancement, depending on the sign and form of the interacting term). Most of this de- 
pletion/enhancement will take place from the background, which in turns translates in 
depleting/enhancing the voids in the large scale structure, which themselves occupy the 
most volume. Therefore, any measurement of the amount of matter in the voids with 
respect to the standard non-coupled model can help to constrain the value of the coupling. 

The matter content in voids has received recent attention (see e.g. Refs. [107, 108] and 
references therein) because it seems that the voids are more empty than expected from 
ACDM model predictions. The work developed in Ref. [108] summarizes very well our 
current knowledge on the occupancy of voids. The local group, mostly our galaxy and 
Andromeda in mass content, is about 20 Mpc away from the nearest cluster of galaxies, 
Virgo, and sits in a relatively low-density region e.g., [109]. Ref. [108] argue that in the 
local volume, at a distance from us of about 10 Mpc, there are too few small galaxies with 
circular velocities Vc below ~ 35 km/s. They estimate a factor 10 discrepancy with the 
ACDM model predictions. Other authors [107] have previously developed models to clear 
the voids of dark matter and thus quench the number of dark halos that could host these 
small galaxies. Of course, another explanation is that the ACDM model is correct and 
some galaxies are simply not visible [110, 111]. Coupled models in which a given amount of 
dark matter per unit volume turns into dark energy, regardless of the local density would 
provide, therefore, a mechanism to clear dark matter (and galaxies) from the voids. On 
the other hand, the void phenomenon disfavors coupled models in which the dark matter 
content per unit volume is enhanced. As an illustration, we use recent void results [108] 
to derive constraints on DMvel and DEvel class II models which will lead to a depletion of 
the dark matter energy density of the universe, in particular, the dark coupled model of 
Eqs. (2.22) characterized by the dimensionless coupling The latter cannot be too large 
and negative to empty the voids completely because, some galaxies with circular velocities 
Vc below 35 km/s are still observed in voids. 

Before getting a precise estimate on the value of ^, one can already impose a preliminary 
lower bound on the coupling parameter by requiring the universe as a whole not to be empty 
of matter in the voids today. Figure I shows the current cold dark matter energy density in 
the universe as a function of the coupling for the dark sector interaction given by Eq. (2.22) 
necessary to fit WMAP 5 year data [1]. The effect of the WMAP parameters uncertainties 
is illustrated by the two lines, chosen to vary along the parameter-degeneracy that most 
affect this figure. This rather simple argument is indicating us that ^ > —0.6, see §3. In 
what follows we will improve this lower limit on ^. To proceed, we will assume that dark 
matter halos associated to galaxies with circular velocities < 35 km/s are depleted by a 
factor of 10 as argued in [108]. Note that at circular velocities of ~ 50 km/s there is no 
discrepancy with a ACDM model [112]. We further argue that such a small galaxies and 
their associated halos in low density regions form recently, when dark energy is important. 
In fact photo-ionization before recombination blows gas away from the halos of small mass. 
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Only those small mass halos that collapse at z < 1 are able to retain their gas and thus 
form stars (e.g. [113] and references therein). With this information we can compute the 
maximum allowed dark matter depletion from the local void and translate it into a bound 
on ^. In fact the velocity function can be transformed into a mass function using the 
spherical collapse (i.e., 35 km/s corresponds to ~ 8 x IO^Mq). At these small masses 
the velocity function and mass functions (both baryonic Ref. [86] and dark matter) are 
well approximated by power laws. We thus modify the ACDM-prcdicted power law for 
the local volume mass function to match the factor of 10 depletion below 35 km/s, and 
still matching the ACDM prediction at 50 km/s, concluding that the total dark matter 
mass in voids can at most be depleted by 20%. This translates into a lower bound for the 
coupling of ^ = —0.2, see Fig. I. This is a sharp lower limit, as it relics on the results 
of Ref. [108], where error bars are way much smaller than the size of the effect. If we 
adopt a more conservative point of view, arguing that the matter is not actually missing 
from the voids, but that star formation in small halos is quenched, then the coupling ^ 
will be limited to the range —0.2 < ^ < 0. This constraint is more restrictive than those 
obtained in Refs. [24,67]) from CMB, LSS and Supernova la datasets. Therefore, future, 
more precise estimations of matter abundance in voids could restrict the phenomenology 
of DEvel coupled models in a stronger way than future CMB experiments as Planck or 
EPIC, sec Ref. [37]. As our calculations suggests, voids offer a promising avenue to the test 
dark coupling. To make our argument more detailed and quantitative, theoretical modeling 
need to be complemented by numerical simulations. 

5. Conclusions 

In this paper we have shown that, although the number of interacting cosmologies that 
have been proposed in the literature is vast, it is always possible to classify all existing 
dark coupling models in two broad families (DEvel and DMvel). Within these families 
we examine two subclasses of models depending on whether the coupling scales like the 
density of dark matter (ex Pdm) or of dark energy (oc pdc)- At the background evolution 
level there are no differences between DEvel models and DMvel models: the background 
evolution depends only on whether the coupling scales like the density of dark matter 
or of dark energy. However at the perturbation level the different classes show different 
phenomenologies which we have explored here. 

GR predicts an exact relation between H{z) and the growth of structure for non- 
interacting dark fluids. Interactions modify this relation, and so can appear like a mod- 
ification of gravity when the growth of structure and expansion history are compared 
precisely. We have shown that DEvel models, where there is no momentum transfer to 
the dark energy rest frame, induce a "fifth force" on the dark matter proportional to its 
peculiar velocity, violating the weak equivalence principle and deviating substantially from 
the uncoupled case in the growth of structure. However, for DMvel models, deviations from 
uncoupled growth arc very small. Finally, DMvel models oc /j^m a-re effectively indistin- 
guishable from minimally coupled dark energy models with time evolving equation of state 
parameter. Therefore, if as a result of a forthcoming experiment, the measured expansion 
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history and growth of structure are in agreement, this observation can be used to constrain 
deviations from GR, as commonly accepted in the hterature. The conversely, however, is 
not necessarily true: a mis-match between expansion and growth could, indicate deviations 
from GR but could also, in principle indicate a dark coupling with GR unchanged. For the 
DMvel coupled models we considered, however, this effect is small. This is summarized in 
Fig. VI. 
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Figure VI: Summary of classification, phenomenology and current constraints on dark coupling 
models. cIl denotes the luminosity distance as obtained from e.g. Supernovae la observations and 
Da denotes the angular diamater distance as measured e.g., by angular BAO. 

We have analyzed how low redshift and near universe probes could be used to constrain 
coupled dark matter-dark energy scenarios. We have considered mis-match between high- 
redshift (Cosmic microwave background-CMB) constraints and low redshift measurements 
of the background quantities such as the Hubble parameter. We have also considered tests 
such as the skewness, which may be promising only for DMvel &; DEvel oc p^e models; 
in these models deviations in the skewness for the dark matter distribution could reach 
the 10% level. It is however not clear whether this effect could be measured or it would 
be degenerate with the dark matter-galaxy bias. For the baryon distribution, (and for 
DMvel and DEvel oc pdm models matching CMB constraints) the deviation barely reaches 
1%. We have also revisited velocity-related probes, such as local bulk flows and redshift 
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space distortions. Present data are not very constraining for DMvcl models, but future 
redshift-space distortions measurements from on-going and proposed galaxy surveys have 
the potential to tighten significantly the allowed coupling window for both DMvel and 
DEvel models. Finally, we have shown that in the context of DMvel & DEvel oc pde models 
with negative coupling, the matter abundance in voids can provide a method to constrain 
interacting models efficiently. This test only depends on the background evolution and can 
be applied to any interacting model depleting or enhancing the dark matter abundance 
homogeneously in time. 

For each of the dark coupling classes we have considered viable and stable models that 
fit CMB observations, and confronted them with present low redshift probes. We have also 
quantified the potential of future data. The current status on the constraints on the viable 
dark coupling models is summarize in Fig. VI. The reported constraints on DEvel oc pdm, 
DEvel oc pdc, DMvcI oc pdm and DMvcI oc come respectively from weak cc|uivalent 
principle violations of galactic scales, voids, redshift-space distortions and voids. 

The combination of precision tests of expansion history (BAO, H{z), supernovae, mat- 
ter content in voids) and growth of structure tests (peculiar velocities and weak lensing) 
together with weak equivalence principle tests on astronomical scales (galaxies and satel- 
lites dynamics) is the key to explore or constrain a possible coupling in the dark sector. 
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A. Background evolution in models with Q — ^Tipde,dm 

For DEvel or DMvcl models with Q = ^T-Lp^ei the Hubble expansion rate function as a 
function of the redshift reads [24] 

H{z) = Ho^n^iii + zr + - (1 + ^)'^- )(i + + ^^£^(1 + , 

(A.l) 

and the redshift dependent equation of state w{z) that one would reconstruct from Eq. (3.1) 
reads [24] 

w 

*W = ; —TTir ' (A.2) 

de 

where = u) + |. If Q = ^V-Pdm then: 



C l^V^' i/o(km/s/Mpc) 






0.1099 


68.29 


-0.1 


0.097 


69.48 


-0.2 


0.083 


70.68 


-0.3 


0.067 


71.90 


-0.4 


0.050 


73.13 


-0.5 


0.032 


74.37 


-0.6 


0.012 


75.63 



Table 1: Current values of the dark matter energy density and the Hubble parameter necessaries 
to fit WMAP 5 year angular diameter distance data as a function of the coupling For numer- 
ical purposes, each time that the coupling is varied, the present day values of the cosmological 
parameters are choosing accordingly to the results of the fit illustrated here. 



H{z) = H, 



\ 



(A.3) 

where = — |.The redshift dependemt equation of state w{z) that one would reconstruct 
from Eq. (3.1), reads, at small redshifts [24] 



(0) 
dm , 



(0) 
de 



(A.4) 



B. Cosmological parameters according to WMAP 5 constraint 

Assuming a flat universe and perfect measurements of ^^^^/i^, O.^^h'^, and the angular 
diameter distance to the last scattering surface from CMB observations, the amplitude of ^ 
is degenerate with the physical energy density in dark matter today ^^dm^^- Table B illus- 
trates for a selected range of couplings this effect on some of the cosmological parameters 
that have been used in the numerical calculations presented here. 
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